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(57) Abstract: 

PROBLEM TO BE SOLVED: To provide a parallel matrix 
processing method suitable for a shared memory type 
scalar computer. 

SOLUTION: A matrix to be LU decomposed is divided into 
a block D at a diagonal part and a row-direction block 
below the D, for example, L1 to L3. D+L1, D+L2 and D 
+L3 are allocated to three processors, respectively, to 
perform operation in parallel. The block U is updated by 
the method of LU decomposition, and further C1 to C3 are 
updated using L1 to L3 and U. The processing is 
conducted for the inside of the block sequentially 
decreased in size, thereby lastly leaving the part 
corresponding to the D. The D is LU decomposed by one 
processor to end LU decomposition for the whole matrix. 
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fr?ij ft /hff « ?P^«i:»lt4?Pv*fcXfv?i, 

KBJi:*-r ^ ?T*» S>*ifc<M¥M ^ p *v * ttBJCBX ft 
Sot, BWTW*I«r-r«E«x^v^i:. 

i^P^ASttttLfc. 1MttBK*ByWBttBB 

tai*a3] «f«a«»fl)^pt7^ya-^*t«i 
jyKD^yBati-fcivc* f->««*is*te 

s^o. . ._. _ _ ;_ . .... 

[M*B4] «TELU»*I4. mffi*T*J(D*«*&W» 

tt<tr;u^ux-arcj:orii«ff«ia>HirtfrLv Mis 

fr*tJOrt. S&(zSgrL»of-6P#ft* l^P-b-; 
oivcLU#»ftKT**C*ft»«fc*4MM2lz 
[MM 5] Kff^J/TMH*. fTMttaUX*— #«fc£ 

«LB-3fcffl#ft* 1 0(D^P-tz^-^^vzL-;U*T?LU 

6 sticky. KffM±f*icoivcLU#jRft3E 
7-r « c <t swat -r 2 

[MM 7] MIEHMx^v^irei^r. 
ffM^n^^(D = flff9Jffl»SI(rlsa»(Z)^P-b^-9-^ 
vzl— ;KD»<Z)2ft(D»(D^P^^lz»aL. R»«S 
titz = fti=r&mft(D7n ^ ft 2 ofoH^tr> 
*?Dt^yi- JKDP- ftJUy^yiftfCftlAL. 
/1IS^^P^^^^V3.-;UIZ^1±^Z<t$^«i:-r 
&Si3fc3g 5 IZEtOffiftHftc 

MMOMMftlBn] 



[oooi] 

[*W©IT4ft«»»] **Bli* **>=EUBX* 

[0002] 

[&*<D#*Pj] B3t1«*B*ttt-JWil=J:-DTB<Ji 
Sit 1 **B*ftfrJUB*U ZflMTflHc-al* 

[00 03] Bai**aai4. <&&ftair 

Ctr. LU^S? (LU Factorization) &1*5:&&I=c*:*l 
I*. fifeft^a^fr^JSJtHftff^J (upper-triangular m 
atrix ) irTHftfT^U (lower-triangular matrix ) iz 
1&MTZ>^t\~£-z>X. 3Ut 1 *»B*0>»ft*tf>4. 

r«ct^iiiti**aaoBt»*fcw=BBtt«B 

(Cholesky Factorization) tl*5fr*iJ<D 

Stfr5iJ(Dilit1**Baa)«^fcBLT. h^UStM 

LU#»ft^— xfcft^HtftfroTl**. M^{7 h;uS 
fl*^*bfc?P.7^SLU»* (1) LfcB, »«L 
b>**tefc.?SL^^.*JE!RL-C (2) 

jE»/hffMtHB o) r**M«iyiB-r. 

[O0 0 4] fit*, (i) (DfliSli. -oo>^Pir^* 

-cB*»icfr-3-ci*fc. M5ija*$*^4fc*!>ic^p^ 

gft^f) B«©lt«M'J*Sl*»i:LTl*fc. C<D$S£ 
(2) attf (3) <7)SU#<dMSt^*Ii 2BB0)fT9«BJK 

[0 0 0 5] «*=iXh<D**l* (3) CDftmT% 
^USxa^afe^JM-JMi (SMP) T?te, I* 

*[ooo6] rftfe^. (3) otntii* fr*"J«T?fc 
4o Ba&«/h*t^H«-r*fTJ«oBB (**"Jl=»B* 
*TCl*4> ftP-KLT. H«re*ftX h7-T4tl^5 
UlcT^-tzX-f £fc#><DaX hA<(TWa)EgrftfT5 
B*l=tt^T**< fcyttBft?l*tt1MPl*. 
[0007] zotztb. 7p?9«ft*«<r«&x*< 
^P-y^«ft*'#<f *tt. :o?d^ol 
u»BOax htf**<*yjWfc»««***. 

b ) 5Efiia»(TM(DjEtt i *»ga©)Ba 

!EB»»f7»©jllt1 **Ba«)BaicBLT. TEft 
f7«»»CD^ft«B LTa UX4~»»ftff 3 t 
»Ky j E , jatSft«ltHlC«-CI*/hff5II^Pv^ftc y c I 
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(DfLffiz&mi-LTMi^T^tzo *fTW coast 1 

SMPtlfc ±C (3) a>S«rt?atu4frM»-e^n^ 

*«#*#i*;&jMt«tfai*fc*K ?p 

[0 008] 

[0 0 0 9] *%ffiO)RHl*. **y^'JIX*7JtI 

acia l fettnfr* Asxa eatt-r « - 1 * * * . 

[0 0 1 o] 

Wf&ztnvtitx'rvzrt. m^'nmz? unvote* 

;uA<stwi=. *ft-F:Mrr a ?p v o zmn? Z>Zt\z 

[OOI 1] *&W\Z&tl\t % ttSlD^Pt^^Sa 
[0 0 12] 

[0013] **y^uffix*5Stww-«w*. 

ZfU-iZVV"\ 0-1. 10-2. • • • 1 0-n^2* 
^t^va^^'J 13-1, 13-2. ■ • ■ 1 3-n 

0-1. 10-2. • ■ ■ 1 0-nli. ^CDflgflfc&l^ 
14. 2**^-r>a> : E , J 13-1. 13-2. - • • 

1 3- n&yyai>yy-m\zi$i*-* visa. y^EU^E 
*^P-b^*1 0- 1 . 10-2. • 

• • 1 O-nC^dr^-DtlN^^^U^va-^l 1 
-1.11-2. ■ • ■ 1 1 - n 14. ffiSfg'&IB 1 2 £ 
^LT^P-fcr^-y-10-1. 10-2. • • • 1 0 - n 
3b<T^*XWte<k*oT5,> ^.Q-fe-v-tM 0- 1 . 10 
-2. • • • 1 0-ntff-: 5»«&S£^5*^l^li. ^ 
"T. ^^'JivaHU 1-1.11-2. •••11 



•r>a>t'J13-1, 13-2. • • -1 3- n left 

[ooi4] «W<»^«fc, i;^tr>a^'j^ 

U ^ya — )\s 11-1.11-2. • • ■ 1 1 - r\ C0 

•b ft #0)5*— * £ 2 X.** v *>=l* * g iz 

oT#T. ^P-b^-9"A<«iS£ff5. 

ft^P-fe ^-9-3b<fflaLf=tt(7)T— * ^ 'J^va — ;U 

^ya-JU:f-^ML^^^t:, te<D :7 P -fe +f 
'J^va-Wbf-^ £ BM*ii*#l*«fc 5 1^ 

5jMM<&*>o ccd^5i^ ^p-tr^-g-^y^g^^i 

^0>ffi8(DBM££ftC££/<y7l^»l (Barrier S 
ynchronization) 2rlfcft <t 5 o 

[ooi 5] H2&i;a3ii. *mw<nmmmiz®.^ 

L U»«fl)15iJfflla)llli:$BMtSBX'*Sc 12 

fr» ^ p ^ (d « sa* ck y * # < u r l u ft a* 

L 1 ~L 3CD35^(4. ?D-^S*'h*<LT, 1 -0(D 
Zfnty* (PE. fcftlM*. XU^K) -CMLtl^ 

tztf. **SB»ttl = eL^T(*. C(D:7p 

-Tfti^lC. L 1 -L 3£-t*L-?tlSlJ<7)XU^ Kic«y 

2*§tt. atOTicjaasi**. fete, cc-eii. xu^k 
[0016] uax*^atMttK«a>4^p-b 

5 - 1 ttffi £ ? i * u -r ± m^-e fe ^ . 

[0 0 17] f-^t0*t4HaS*»P ET*«<li 
^. ftp ETff- ^^^bLT±i*4: LTIi^P^^ 
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[0 0 18] ZOtzM^ m3\Zijk£tlZ>&ol~* 

45i)l:I «t 4 1 f MS ? □ ^ ^ »ft 9» D li, 

fit?? h£S*ofc&* fr** hJU(DAtil§;L£fr5i§£- 
14. ivfjh,J&*o>:?n-fe 2 * 'r>i^'J ±1= 

[0019] b : ^D?£flL k : *^a-tz 

ft* *T>tft>*». H2l:fel+4L1-L30?P7^0§ 
fKDftCDft) £Lt-<b£\ bx (b + k) xs~8Mb. 

[0 0 2 0] *LT.. fl?*« (l^^vi^U) 

*£WJBLfcLU5MR*fr5. ««-r^#ff»lafi(< 
^ 'J {ft&\tf*i£<) . **Llrtfc 

^<7(D^p ^^«3b</h* <-65fc #14* :?p 
[00 2 1 ] Mlc *^P-b'V^ («l>7 K) T'CDfE 

[0022] 114(4. *sisfi^»a)Lu*MBa)ffiaa)3SE 
.fte^rwifc^p— htfe^o *-*\ ^f7?s 

i) ^i^l:^t-t4o f LT» X^^^S3lcfc 
1>T\ ftXU7 Ktftftf? 

K-Cr. ±E?P79Dfc?D79L i fc*LU#IHt 

[0023] ^t-.^^s Aizfci^ii. ffla*<»*?y^ 

f**>y-efci**fti;m zf7^s smfci^T. 

KT?365>JIc^p 7^LL (B2#IR) JS^t, :7p 
•^Ui (B2M) ^ LU ^S©7JHJ U XAII^o 
tlit^o fit, X-TV^&6lCj3L^r. «l/7 
K-CatWJ-. ?P79Ci (B2#SR) S^7H.i 



t^P7^U (U i t>l£tz*><D : B2M) 0) 

mx-'g.Sit&o fit. Xfv^S2i:RU, XO^p 
(Ci^!»447P7^) B»©»atLU» 

«-t*. tut, *Kic/j\a<ft«*»«SDv^*<* 
8biz^p^^Di=»B-r«»#a>w=<Ey. ioc^xu 

7KtLU»S^7f5t, fr5IJ±*lCO^TLU» 

[0024] 15-11011 *SH6^JB0)LU»»(D 
3Sa**yiMBlcl!lW-f «B"tf*«. -:t*lt 2 04 
8 x 204 8fl)W)$4Xb7 K"C L U»»"T 

[0 0 2 5] ff, 2 0 4 8 x 204 8(DWJ^^P'> 
9«2 5 6T'LU#«t4t(DttSo ft?P7^(Dg 
#1+14. B5l^1"t fcyCfc*. ibt, 40CDCP 
U WU7F) l?«at*fr-r*«*. SXb7Ktl 
tt«;1l«t ((256 + 448) X256:8MB (L2 

K^-r<fc5C^ &XI/7 Kfl>#fHW-D 1 +L 1 . D1 
+ L2. D1+L3. D1+L4S3t-t4o 

[0 0 2 6] £*S. ^P 5. ttTa>*5l=' 

LT»4o RSH<&*:#* (ffMO)*** : ffWCD« 
ft) £n. Xb7K<DSS#THREADi:Lt, 

[0 0 2 7] 

[ftl] 

^(n*/# THREAD) x^==nb 

[0 0 2 8] tfcl^X. 
nbi512ft&, ^D7^fi=512' 
nb2256&&« ?D7^B=256 
nb|128$b, ?P79«=128 
**tlcl*k ?P7*i=6 4 

[0 0 2 9] t6fc%, LU»»(DaXHt 2n3/ 

t % h^#THREADt^lHbl, ^0)1 

%*<iilw17l/^ K"CfT3^P^**«* £fc£><J:5 

[0 0 3 0] CwT*. 3a»SHjlt*fcto. SfcJMtLttl* 
S^<DLU»«(D7;Ud'yXA$B7lC^-r o B 7 left 
LT = D1+L1+L2 + L3 + L4^J 
LU^t^T^il'jXAS^Lt^^o LTI4. 20 
48 X 2 560^D7^ (h^otl^e 

[0031] *-*\ (i ) oajtficfc^T. h£ 

aS"T5o iblks l*LTto«tffcy. ^(D«*. 2 56 
"CftS. *fc* Iengl4. LTC0«*-C&y, ^CDi§-&. 
204 8T'^p j j 1-14. tf#? hOffS-TSfr** 

[0032] fit, (2) o>»»icJ3i\t. aasaa 



(5) 



ftffl2 O 0 2-1 63246 



[0 0 3 3] C(D ( 1 ) ~ (3) £ i 7b<1 ~ i b I ks {eft 

/=oT«y«L3lJcr4o ::t\ LT<7>ft£T*fc£. 

Ieng*<i^<h*£<£££. Cti&OfflSI* L 2*-^ 

T£3l£®-*To f:t\ E 5<7)cfc die tF— £ 

L. L2***;/vi>^yi-^— ££RI#Lfc££4l&3 

[0 0 3 4] H8lzj3L^TIi, LT i I*. P — * 

;u«L pivot (4). GPIVOT, ROW (iblks ) It. ^fflHElc 
te^£ft£^--*T3b6o *"*\ (4) lcfcl>T\ ftx 

K-etr*^ h£i*£o tit, (5) -e«*tt£E 

5iJpivot(4)(DXU^ K*#0)S*lc»lrt-rSo (5) 0) 
A»JrB#i£S&oT. (6) Izfcivc. fi*tftf 
^ htl»OZU7 K»#SGPIV0TlC»«-r4o fit, 
(6) (D&T*. Stf/<'J7H«IJ»«. 2W=. (7) l= 

(8) (Z*3LVCI4* GPIVOT*<0 0>££I4. **fcf#*h 
I4D i <Dto\z&z>i)\ Ah»^5t-fcy, p-*;uiH 
lcA*t^*-£o GPIVOT*<#T H READI^Ll^#, 

■rtr*>%. GPivowocfcy**L*Jt*i=i*, 

h;u£ A*ilf*-£o tit, (9) . (1 0) |=Jsi* 
T\ LU#»<&fca&(Z>3m£fr5o tLt, ±E (4) 

- do) SffiS-r&^n^a^TtBMlzo^TfT 
5o "Tfe^^. 1 -iblks £T!CD i (zo^T«131£l*y 

[O 0 3 5] CC-C. LTi0LU^it-C/\'J7 

1«-LL-1U1* U2<-LL"1U2. U3<-LL^U 
3. U 4<-L L-1 U 4 K-CM^lJlcSta-T^o 
CCDf+glCD&T;. D1, LI, L2, L 3 „ L4$D — 
*;U**&fr*]Alca fcf— L. /^JTHMfR 
£«> Ml-. CK-C1-L1XU, C2<-C2-L2x 
IK C3<-C3-L3xU % C4<-C4-L4xUfl 

?|JI^S< K"tffrl\ *«i:/^J7H«li»S. 

[0 0 3 6] S9I4. ±Effl«l-J:o"C* ISBOBI 

H9lr5?t"«k3lC % DK L1, L2. L3. L4CDZ? 



C2 5 6 x 2 5 6fl)?P w ;^tfS4. C0)»»Ji-1O© 

[0 0 3 7] fcfc. ±Eftai?l*. LT i jLU^gt 

6fc*IW»tt<tLU#»ft«fflLT^*. C i<- 

C i - L i x U <Z)3Mtl4* ^ ^ ^>zl ^t'J JiCD-tt— £ 

[0 0 3 8] HI 014. I»«LU»(g7^»jXA.$ 

wwttLutfwoTiu^uXAtt-y- 

^U-f >LU<tLT4^btll)o L UCDSS££$fcl4. 
BQflDT^yXA-CBt**:^ LT i KT* 
D 1 + L i Sfttt) , k (LTi a>1*5cB0)** 
S) „ iblks (^P^^«) COtelz. LU»«**6tf>4 
tem£^*Tist. LU»»*?75(I"efc5riwid-efc5a 
[0 0 3 9] ^><7>fit0-e. nwid*<8. 

£ P Y E S<Z)»*ICli, LT i (ist : k. ist : ist 
+nwid-1) *a6MlcLU»j»*r*o ist : k 

irs^aisi*. ft»*<ist ^5,k*-co)LT i £^-rm 

tt-C. ist : ist +nwid- 1 <hL>5a>l4. fcHtfistft* 
&ist +nwid- 1 *t»<DLT i $^-T*i*-Cfc* 0 KIT 

[004 0] LT i (D.LU#»lCfel\TI*. ±E < 4 ) 
-(10) MIiit5 fi fc*£U fTa>A*x»*.»# 
14. ft* iblks T*. LTi (u 1 : iblks ) ^Atif 
*fc % ±E*JWTA<NO0)*«lcli. LU#»Jft 
5»Hnwid£ 2 -CMofcttO) L U#j§?£<D L U^P^if^ 

T>&Wlfft-?o ZCDJl—^yi*. LT i (is 
t : ist +nwid/2—1^ ist 4-nwid/2 : ist +nwi 
d) £MST-r*o MIC. LTi (ist : ist +nwid/2 
-1. ist : ist +nwid/2-1) OTHflffJULLS 
Silt, LL"1$C i lcS^e>*MfTS«r-r*o * 
MMi:l^5;^— ^^^Dft/a-To Z<DiU— J 

LT i (ist +nwid/2 : k. ist +nwid/2 : ist + 
nwid) = LT i (ist +nwid^2 : k. ist +nwid/ 

2 : ist -hnwid) -LTi (ist +nwid/2 : k. ist 
: ist -hnwid/2— 1 ) xLTi (ist : ist +nwid 
/2-1, ist +nwid/2 : ist +nwid) 

Jllt^o fit, /<'J7B«**y. LU 

b) iEfllWWfTWCDaEa i *^a*(7)»a 

Si i -mi 3 14. M»»f7»0)»#i:al/X*-» 

[0041] nff^jcD^iL i fr?"j$ 

D. LI. L2. L 3IC»«ILT. *XU^ Kir*fftfr 
^JDi:M«r-r«>OT^P^^fiP»L 1 . L2. L3SM 
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Stlczi tf— LT (HI 2#H) . &±L\z&5\\\zLT&)7 

*MLT. 'hTHflWJ (C 1 ~C6j^bfc£) £M$r 

Z>tzib\z< HS»-r*^#T = ft»5iJSXL/'v K»£#T 
fcLfci*. 2 x#T|zHi:^P^^Mlw»l+4 CTtt 
Cl£C6. C2<tC5, C3tC4ji^ 
■&-C»a*R5J:ai=-r«) o «U*; KI4 i *SAlf 
2x#T+l-i SKD^a^^SSSff 

[oo4 2] aa*5ir^*ftW3a<iE«»iMT«-cj54 

CD"C. H2iDUlZ$*J6-r4»»l±. L1+L2+L3^ 

^lc(4. lT as#i4. 3W£ir5i&KjWR<. Li. l 

2. L 3 ^-^^-rtvfeMLrzi t°--r*t!^J:L>o 
[004 3] 11 3li, lsX*—&1ltO>8!M(0 

T. fTM<7>— S£CD^P-y^S3 U^*-»*L, 3 0 
<D»#£3 1 ic^e— "T ^LT. 3 1 <DT<Dj£fg-T? 

(2) irfcl^T* 2fil: LT 3 2<D« 

#£3 3 1:: =a tf— L. 3 3COTCO,#^T*S*tt/ca5^^ 
»«»»*Mt^T. WfrT*. *'lt, (3) \zmSu 
V. 3 4<Dffl#£3 Siratf— 3 5 0TO««t* 

»«»»sffli^rM«-r«« sic. 

(4) icfc^T. 3.6JC»*4x4»»* 3 7 Jf— 
L. 3 7a>T0)»#t*Hi»»Sffl^TE3R"r«. s 
(~. (5) IZfcnT. 38^39|:n tf — L. 3 9 0)T 

(6) iZfcUT. 40$41|C3tf-L, 

l N *C. 4 1 ©T©9^$I8rtl)« (7) lcfct> 

T. 4 20)»»S4 3lZ3tf— L. *HB»»*ffll^T4 

[0 0 4 4] HI 4 -Hi 614. SE^P UX+-»*0) 
14. lttW£«Kfcr«fc*>4XU? K$«ot*«fi 
[0 0 4 5] KI3H 1 4 (Z^^tl^) D x 

&tf l i $ian«««LT » icatf— r*. *ct\ 

.LTiSLDLT L D LT #g!4IS»ttfc* 

a-Cfr5. fit, DL i ^iaTOftg-e<!£ft*iJI::p 

if-**. 

[0046] DL i «-D x L i T . ::t\ DliDxO 
«o fit, Dx (7l/7Klf^b) te<7) L i I4££ij 

ir. i>ta>nmzzi tf— fit, ::t-, 
/«U7H»*B l J . Hi 5ie^£;tx£cfc5ie. CUC 



8. C2<hC7. C3«tC6. C 4 d: C 5 ?^7<i: L 

C 1 «-C 1 — L 1 1xDL1 1 
C 8 «- C 8 - L 4 2 x D L 4 2 
• X U'> K2-CI4, 
C 2 «- C 2 — L 1 2 x D L 1 2 
C 7«-C 7-L41XDL41 

K3TM4. 
C 3 4- C 3 — L 2 1 x D L 2 1 
C 6<-C 6 - L 3 2 x D L 3 2 
•XU^K4TM*. 
C4<-C4-L22XDL22 
C 5<-C 5-L31XDL31 

[00 4 7] fit, zc*t?a>3SK*<lSfe-3fc«jSn? 

[0 0 4 8] ±f£. ft*U* Klfff5C i OIMIS 

cfcy-»ttft»K-caL-<*fc. xi/7K»s#thre 

ADt LT. L£2 x#THRE ADffllC^tlL. C<D 
T = ft»»tHl^< 2 x#THRE AD^iC^flJ-r^o 
fit, ±tT*b#THREADlO)^7jfty, - 

(D^7t c(D»»j Lfc»#*3E»r-r « fc 5 naictt 

[0 0 4 9] H1 6I4. ldlT »»©W»tt7^U=fU 
XA^titfcl) 0 L D LT- ^©Vw'ixittJ' 

Bi££gfcl4. M(0LU»(Bfl)«4i:H«tfc6o 

[0 0 5 0] £"f\ nwid*<8(D**iri4. (2 0) (DSP 
#T*. It^LDLT ^SStr5o fits LT i (ist 
+ 8 : k. ist : ist +7) £H£rt*£o CKDfc*. L 
T i (ist : ist +7. ist : ist +7) CDJLH^ISJW 
COLT tfA-oXl^&OX'. (DLT ) ~l£:&fr&ftM? 

[005 1] (20) <7>S*J[)CD I F^tT*. riwid*<8"Cft 
l^M»r S4xfc»«icl4. nwid/2 £6r*=£nwid£ it 
^;u-f>LDL$u?i/t±iL, HfT-T^o ::t% l 

T i (ist : ist +nwid/'2- 1 . ist +nwid/2 : is 
t +nwid- 1 ) ICD LT £=j tf— DI4. L T i 
(ist : ist + nwid/2-1. ist : ist + nwid/2- 
1) ©SMIWircfcy. LI4, LT i (ist +nwid/ 
2 : ist +nwid- 1 . ist : ist +nwid/"2- 1 ) Tfc 

y. :(DLjgit^ c 

[0052] fit, LT i (ist +nwid/2 : k. is 
t H-nwid/2 : ist +mrid- 1) £M8r^£<> -Tft*? 

LT i (ist +nwid/2 : k. ist +nwidX2 : ist + 
nwid— 1 ) = L T i (ist +nwid/ 2 : k . ist +nwid 
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/Z : ist +nwid-1) - LT i (ist +nwid/2 : 
k , ist : ist + nwid — 1 ) x L T i (ist : ist + nw 
id/ 2-1* ist + nwid/2 : ist + nwid- 1 ) 

[O0 5 3] LDLT ^SOf^-f >LDL 

znffiifti~wit&To fit, mmtfrnfr-Dtz*. *yy 
jis-^^zfott&o tet$s *§%*m^mm\±. ±ek 

r/up'JXAfc it4^.bti§a)-e, zcdtjup uxa 

ffl«fc L-Cfl6ffl-r4»-&lc|* % ROMfcfc'lc^P^A 

ttfc LTffiffl-r*»*i-i** *xna>xifcjiMBa>rju=f 

yXAIi. CD-ROM^<DBr»E»tti*-V>. /\— K-r 
-f x*tt£a>E»MM*lc3fay^Afc LTEBLTfc 

[0054] c<Dj:3<i:a*. ^mmmmmoT )i, 

=TUXA«|U|-r ft^ay^AI*. ^*E*«<*fc£S 
Sot, a.—- tPlcE*A<RTte-efc4 0 
<•**«) 

• l u#M<D;£8fcii> i) % 2) % a^ux^S 

<7)^ffili2) 

1 ) P. AMESTOY, M.DAYDE, and I. DUFF, "Use of c 
omputational kernel sin the solution of full 
and sparse linear equations" , M. COSNARD, Y. R 
OBERT, Q.QU1NT0N. and M. RAYNAL, PARALLEL&D I STR 
I BUTED ALGORITHMS, North-Holland, 1989, pp. 13- 
19 

2 ) 6. H. Go I ub, C. F. van Loan, "Matr i x Computat 
ions" , second edition, The Johns Hopkins Un 
iversity Press, 1989 

• B*B(D*i»-ei4. JaT<Dt<Z)lCLU#«;&l/LDL 

• rRffigWrj ft IEMK ^£tiiKfc:£*± 

$q ? a * |z» It « y u v >7 <bx t- ? ^ fc . M'Mt*'J ? 

fc fiS»(D^Pt7«vaHKDD-i3^ 



^A£te*ALtr. nmMRMBtj pTfig^E^Kfto 
[00 5 5] (W-E 2 ) ^^T?lJ>H» (4. fir Jfl<0 L U-»» 

■c fc 4 c fc £ fc r *<*e i i:e8oe*«*. 

fc**f*E2lcE«©E*«tt. 
[0 0 5 6] WE 4) WELU»«I*. WrE'fTMCT)^ 

frlV «rEff3W0)rt. **lcE*L»-Dfc»#£* n 
0>^P-te^-9-^5?a— JU-CLUam-riCfcfCcfctj . K 

^j^IwOint Lussj^Tt^c fc -r* 

«'E2lcE«0)E»JKf*. 

[0 0 5 7] (ftE5) mft&mm*. frJljOP bX* 
--#Sfc4lMMM£p l/ , X+-»*Tffc4 Z fc £«fc»fc 
"T SttE 1 lcE*0)E»«tt:. 

(f*E6) WfiaUX*— #*U&4lM«»aU;**— 
lWE*TW©^«^6W»»ttTyu=f'JXAlc* 

PtitilictU, Kfffl^ttlcoivcLUtflKSSST 
f « C t * W« 4 "T * #E 2 lrEKOE»«<*o 
[0 0 5 8] «tE 7) KrEIftX-T'^^fc^-C. E 
gr-T # /hfT*J ^ p * O = ft ff5"J«» * ^TES^IO ^ 
P-b -y^^Ei/ZL— ;UG>ft(Z> 2fS(D»(D^P ^ iC^flJ 
L. BE»**4xfcHft(TWM»a>^p^*2oro« 

L. ^P-fe S/i-JbJCfiW* c 

fc * »« fc "T * ttE 5 lzEKcDES«(*c 
[0 0 5 9] (WE 8) «»0>^P-b^^^i-;U$ 

r , ff « t /M¥*j ^ p v * s ^ p v >? {tTs^f v y 
fc. m^?fm7av<?oft. ttnyo-yotttnyny 

tX'tZ t ^ P * fc ?Pt7«va 

-;kdp-*;u>^ ( J «*lz»*-r**lftXT^^fc. 

4^p^**aijnr6zfcic«fey. a»a>^p*^** 
yi-jbT?, wft^p^^^jcfticajt-rsatjr^-r^ 
^fc. ttajrx^^^wfetifc/hffM^n^^oatji 

«ct**tm it jiff 

[0 0 6 0] (f*E9) MO^Pt^T^va-^ 

fi- J>J ? P v <? Z> y p v 0 fc . 1S/J^^TJ^I ^ p 
^^©W. Wft^P^^fcWfl^P'y^-Cttlvhffjij^ 
p fc Si$a»i7)^P-b ^"9-^ va-JKOP- *;U* 
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& tufc'MrM ? p v * Jftot, BHtW * JE 

A1S1. 
[006 1] 

[H2] *«wo*it»tt^«ofcLu#»a)SkWffla 

[13] **WO>j|Jfi^Blz«EofcLU»»a)3fe5ilffia 

(-to 2) r*fe^o 

[15] *Sllfi^ffl(DLU»»(D^a*J:yP«lcKW 

[16] *as*»«a>LU#«(a)*ft*«fcy»«Bi=lftW 

[17] LU»fKD*ft*cty»«l^ttM 



-f£l (*0)3) T*fe£o 
[18] ***»S0)LU *y»«KBWS 

[19] *3EJ6^ffia)LU»»a)^ra^*yBffll-HW 
[li o] **tfi»B<7>LU»»©*a*«ky»ttlcK 
[111] jEffl»»fT«0)»*l=aux*-#i(*fr5 
[112] jEttWttffWO)**!^ UX^r-»SJft5 

[113] iEmttWnm<nm^{zz3[s**-i*Mir?Tv 

[11 4 ] X»3 i/^*-5M»©7^^ 'JXAf* 
[11 5] »»a L/^^r-^»CDT;U=f'jXA* cfcUB 
[116] SEJUa U^^-^S©7^ 'J XASi: y B 
[»*<DKW] 



1 0—1 — 1 0 — n 
1 1 - 1 ~ 1 1 — n 
1 2 *§5*S^8i 
1 3—1—1 3 — n 



[11] 



[12] 



IH H-2 



10- 1 



10-2 



13-n 



tO-n 




12 









Nil uu 

u/V 


1 u : 
ui ] U2 ; U3 




LI 


Cl 




L2 


C2 




L3 


cs 
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[03] 



[0 4] 



2^fi*Hfemm\'&*rzu)'jMfr'> Anna. 



D 




D 




0 


L1 




L2 




L3 ' 




[05] 



(7 START 

jzrzi 



3: 




c 



3 



«48 

V 
/■ 

44B 

V 

/ 



Dl 



L1 



L2 



L3 



L4 



U1 U2 U3 U4 



CI 



C2 



C3 



C4 



[06] 



2K,K*1 XUv>'*2 XWr^tS XUvF"*4 



01 




Dl 




Dl 




D1 


Lf 




L2 


LS 




L4 



00) 
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[17] 



[B8] 



DO i=1.fbtk» 
TMP—OlO DO:i}=0 

DO jexu lens 
XFCABS LTQ. J». 07 
TW=ABS<LT0. 0) 

ENDIF 
ENDDO 



IF(iJ. ST. I> THEN 

DO HI, folks 

TMPX-LTU 10 

LTa k>«UT(5. k) 

LTCfi,k)*TMPX 
ENDDO 
END IP 

DO fcBf-M. Iblks 

_LT0. k)«LTCI. I) 



TMP/THEN 



] 



(1> 



do tel. ««■ 

TM ffc=O,0 DO J=© 
DO 1=1. 
TIAftS LTi{j. 0). GT 
TUP-A&BU-T1U. 0) 

EN OCT 



TMfOTMEN 



] 
] 



tF(«rrvREAB. ea. 1> 

DO * 

Flpraot&c}. OT, t»- AND. P3VOTUx>. OT- MnJUWUIi-U 



end r 



in«TTHR£AD. ea. cp*ot)tmbji 

IFCI. GT, DTHEM 
DO br-1. M«* 

ENODO 

END r 



C7> 



C2) 



XFH5PJVOT. EO. 03TWEX 

rti. err, othb* 

DO " ~" 



LT-,0. bO=UTiU- *»1 



una,**) 



DD t«— 1+1. Ibfla 
DO Iftnc 

LTD. W«LT(I. k>-LXa i) NLT0. k) 
END DO 
ENDDO 



D?<»T>mEAD. ED. 
DO 

wnu w-Rowtw 



<3) 



uui w-rowcw 



[@9] 



T 

296 

s. 

SM 

V 

364 

r 

394 

r 

364 



01 



LI 



L2 



L3 



L4 



U1 U2 U3 U4 



C 1 



C2 



C3 



C4 



O] 



subroutine UJtLTi. k. folks. W. ««rfd) 
(LTttt*3lU"yK*rD1 +L*fM6. 

ffHnwtd. *q. a)TW.(rPBrjl / >> 
LTlG»t:k. i*t l**+nwW-l)*JI 



L *^f^DAK»X«*ti**!bikrLTia 1 . folk.)****** J 



csfi LUCLJi. k. ibftui. 1st f»*W/2) 

csfl TRB( ) _ 

LTlOsrlsH-rrvwd/2-1. t*t+rnw^2:l«t+fT»fId) 

call MM{ ) 

CLTia«H-mwkl/2.:k. it+nrnixlSZ : Ist+nwW) 
»LTt( hJ. I iww i d^2:k. tat 4 nw k t^&:t»H-fwrtd) 
-LTlGaH-nwkt/2:k. im»:l«t+rmW/2-1) * 



] 



can LLKLTV k. folks, i s l-i m w d ^ nwid/2 
s«d if 
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[i1 2] 



it'll 5W<??'l « Jf^it a S, 





(12) 
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[m 1 5] 



DU2J 




II 6] 



£9ittticK«i4-«ac*03> 



■ufaroutin* LTD4LTL k. R>»<&. 1st nmW) 
lECnwid. EQ, 80 THEN <lfcB6«'» 
DOI«i»t <st+7 
DO*=i+1. lst+7 
LTi(i.i)=LT.aO 

una o=LT»a o/ltiu. 0 

ENDDC 

DO Jy=l+1. tet+7 
DO J*=*x. fafc+7 

LTJ Ox. jy) =LTi {jx. jy) -LTi (K 0 * LTi 0. br) 
ENDDO 
END DO 



c 



J 



LTKLTi(I>ft+8;k. |«t:i«H-7)«Ur 

U1U.TKlat:lst+7. l«t:lBt+7)0-fcHMl=DL T *A^-a^O-C 



D 



cmII LDLCLTi. K. Ibtka. 1st nwldy'2) 

•LTKi«t:l«t+nwklXg-1. ■ a t+rmklX2:iirt+nwiri--n 
r=DL T eatf--r*. CD<*LT]Ctat:bt+nwid/2--l. tat:fart-rnwid/2-l) 

LT)Ci^+nwid/Z:J«t+nwid-1. ) 



'LTiGaH-rwid^k. la*+T»rid/2:i«+m«d-1>«aS 



LTI 

=L 
LTi 1 
ls_ LTi 



LTKiB*+nwkl/2:k, i«H-nwid/2:i«H-rwrid--1) 

LTI (i*t :«t+fiwW/2:k. l«t+m*id/2:f«H-r»wid— 1 )- 
LTi(l«*+m*idX2:k. fc*:iBt+mwd-"l) K 
LTi(wt:i«t+r«wld^2— 1. W+nwW^^lsi+rwrtd— 1 J 

CALL LDL (LTi. k. ibttc*. «t+nwrid/2. nwid/2) 



J 



RETURN 
END 



